Quantum nematic as ground state of a two-dimensional electron gas in a magnetic field 
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We study the ground state of a nematic phase of the two-dimensional electron gas at filling fraction 
v = 1/2 using a variational wavefunction having Jastrow pair-correlations of the form n i< j(z i — zj) 2 
and an elliptical Fermi sea. Using the Fermi hypernetted chain approximation we find that below 
a critical value of the broken symmetry parameter, the nematic phase is energetically favorable 
as compared to the isotropic state for the second excited Landau level. We also find that below 
a critical value of the layer "thickness" parameter A (and in the actual materials) the quantum 
nematic is energetically favorable relative to the stripe ordered Wigner crystal phase. 

PACS numbers: 73.43.-f,73.43.Cd,73.43.Lp 



c3 



i 

-a 
c 

o 
o 



> 

(N 

(N 
O 

o 

-I— > 



i 

C 

O 

o 



During the past two decades, the quantum Hall ef- 
fect (QHE) has been one of the most intriguing research 
topics in condensed matter—. More recently, the mea- 
surements of Lilly et. at 2 and Du et. al.~ reveal a 
strong anisotropic behavior of transport properties of 
electrons for the half-filled Landau level (LL) system un- 
der strong magnetic field and at very low temperature. 
The anisotropy commences at the second excited LL and 
persists up to the sixth excited LL. The sudden exhi- 
bition of large anisotropies of resistivities in clean two- 
dimensional electron gas (2DEG) suggests that there is 
an unknown underlying microscopic origin for this kind 
of spontaneously symmetry breaking. 

These experimental findings prompted several inter- 
esting theoretical proposals which attempt to explain 
the observed anisotropic behavior of the half-filled LL 
system. First of all, these anisotropic transport prop- 
erties are consistent with already predicted stripe and 
bubble charge-density-wave phases found^ by means of 
Hartree-Fock calculations of the 2DEG. However, Frad- 
kin and Kivelson 6 suggested that the anisotropic trans- 
port might be due to a stripe nematic phase of the 2DEG 
in a high magnetic field. This point of view was inves- 
tigated further by Fradkin et alZ where a model for the 
nematic phase in a symmetry breaking field was studied 
using Monte Carlo simulation. The results of the Monte 
Carlo simulation provide a good fit of the experimental 
data of Lilly et al. 2 . This simulation suggests that the 
nematic phase might be a good candidate to explain the 
anisotropic behavior observed in Ref. and Ref . 0. Fur- 
thermore, by deriving a long-wavelength elastic theory of 
the quantum Hall smectic state, Wexler and Dorsey^ have 
estimated the transition temperature from an isotropic 
to nematic phase to be of the order of 200mK. Later, 
Cooper et alZ by applying an in-plane magnetic field 
in 2DEG samples which show the above anisotropy in 
transport, give further support for the possible presence 
of such a quantum nematic phase. 

In the composite fermion theory given by Jain^, the 
fractional quantum Hall effect (FQHE) is interpreted as 
the integer quantum Hall effect of composite Fermions. 
Furthermore, Halperin et alr^ developed a theory of half- 



filled LL system, which is the case of our interest, as a 
compressible Fermi liquid. Rezayi and Read" proposed 
a ground state wave function for the half-filled LL system 
having the Jastrow-Slater form as follows: 

N N 

* (ri, f 2 , ...,r N ) = P J] (zj-Zkfe- 1 ^^ 2 

j<k 

x det 1^(^)1 , (1) 

where tpzifi) are two-dimensional (2D) plane-wave 
states, and P is the LL projector operator. Here Zj = 
Xj + iyj is the complex 2D coordinate of the j electron. 
This wavefunction is a Jastrow correlated Slater deter- 
minant with Jastrow part similar to the Laughlin state. 

Ciftja and Wexler— used the Fermi-hypernetted-chain 
(FHNC) approximation to study a broken rotational 
state of the half-filled LL where the symmetry breaking 
parameter was introduced in the correlation part of the 
wavefunction and the single particle determinant of the 
wavefunction was characterized by a circular Fermi sea. 

The ground state wavefunction for the nematic state 
proposed in Ref. Ql] has the same form as the wavefunc- 
tion given by Eq.[T] however, the single particle momenta 
form an elliptical Fermi sea as opposed to the circular 
Fermi sea. In this paper we will use this wavefunction to 
study the nematic phase. The broken symmetry parame- 
ter in our problem is the ratio a — kx/k^ of the major k% 
and minor hi axis of the elliptic Fermi sea. We will study 
the nematic state of the half-filled LL system using the 
variational approach and we will employ the Fermi hy- 
pernetted chain approximatio n 19 : 20 ' 21 : 22 . We adopt the 
ansatz for the ground state of the nematic state proposed 
by Oganesyan et alJZ as trial wavefunction in our vari- 
ational calculations. Namely, we investigate whether or 
not this state, in which the anisotropy is due to an ellip- 
tical Fermi sea, can be energetically favorable relative to 
the isotropic state and the stripe ordered Wigner crys- 
tal at high LL. Wc find that the nematic phase can be 
stabilized against the isotropic case beyond the second 
excited LL. In addition, we have compared the energy 
of the nematic state to that obtained by a self-consistent 
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Hartree-Fock calculation 2 ^, of the stripe and bubble states 
of the 2DEG. We find that there is a transition between 
the stripe ordered ground state and the nematic phase as 
a function of the parameter A of the Zhang-Das Sarma 
interaction^, i.e., the layer finite thickness. In particu- 
lar, for the case of the materials^ 3 - in question, we find 
that the ground state corresponds to the nematic state. 

The FHNC formalism for Fermi systems was intro- 
duced and developed in Refs. EHHllil A significant 
advantage of FHNC over the variational Monte Carlo is 
that FHNC does not suffer from finite-size effects. In the 
current problem, where we need to estimate small energy 
differences, the role of finite-size effects may be signif- 
icant. The main idea of this method is to expand the 
pair distribution function in powers of the density of the 
system and it works well for a low density system. The 
variational Monte Carlo becomes advantageous at rela- 
tively high density, as in the case of liquid 3 He, where 
the role of the elementary diagrams and other non-trivial 
many-body correlations need to be included.— 

In order to apply the FHNC formalism, first, we as- 
sume that the unprojected wavefunction is good approx- 
imation to the real wavefunction for the LLL. Moreover, 
it is known 1 ^ that the projection operator almost elimi- 
nates the high LL components of the wavefunction. The 
potential energy of the high LL can be expressed^ 2 - via 
the pair distribution function of the LLL using the single 
mode approximation discussed in Ref. HH namely: 



V L 



£ l[9(r)-l]V$(r)d\ 



(2) 



where the effective potential V^J (r) at Landau level L is 
the convolution of the effective Zhang-Das Sarma (ZDS) 
interaction^, V(r) = e 2 /e\/r 2 + A 2 , with the L-orderLa- 
guerre polynomials; namely, it is the Fourier transform of 
V$(q) = (27re 2 )/e 9 exp(-Ag)[L L ( g 2 /2)] 2 . In the above 
formula A is a length scale which characterizes the con- 
finement of the electron wavefunction in the direction 
perpendicular to the heterojunctionp^ 

Our calculation proceeds as follows: First, we will cal- 
culate the pair distribution functions for isotropic and 
nematic states with different values of the anisotropic 
parameter for the LLL using the FHNC approximation. 
Second, the interaction energies will be calculated via the 
pair distribution functions by using the single-LL approx- 
imation, i.e., via Eq. [51 Next, the kinetic energy is eval- 
uated for the isotropic and different nematic states. The 
energy values of the isotropic state are compared with 
anisotropic states for the lowest, first and second excited 
LL to find out if the nematic state becomes energetically 
favorable. Finally, we will carry out a self-consistent 
Hartree-Fock calculatio n 23 ! 24 for the more general case 
where A can be non-zero and we will compare the en- 
ergy of the nematic, isotropic and stripe ordered Wigner 
crystal. 

In the FHNC technique, each term in the expansion is 
represented as graphical diagrams with well-defined topo- 



logical rules. There are nodal, composite and elemen- 
tary diagrams. In the FHNC/0 approximation, which 
neglects the elementary diagrams, the pair distribution 
function is obtained by solving the FHNC integral equa- 
tions given in Ref. 21 for the case of polarized liquid 
3 He. These equations require as input, (a) the pair cor- 
relation (or Jastrow) factor, which in our case is given as 
/ 2 (r) = exp(ufr)) with u(r) = 4Zn(r), and (b) the statis- 
tical exchange factor l(r), which for a 2D Slater determi- 
nant is given by: l(r) — C(kpr), where C(x) = 2J\{x)/x, 
J\{r) is the first order Bessel function, and kp is the 
Fermi momentum of the isotropic state. Alternatively for 
the anisotropic state having an elliptic Fermi surface with 
major and minor axes k\ and fe, we find l(r) — C(X), 
with X — \J (k±x) 2 + (fey) 2 where x and y are the coor- 
dinates of r. 

Since the pseudo-potential u(r) has a long-range loga- 
rithmic form, we follow the standard procedure used in 
Ref. [lCj and Ref. [l| to separate it into a short-range and a 
long-range part. This leads to a new set of FHNC equa- 
tions for the short range nodal and composite functions 
which are solved iteratively using a combination of mo- 
mentum space and real space approach^ 



•i 0.05 ■ 




>-«L=0 

S-B L=l 

mL=2 




FIG. 1: Top: The potential energy difference AVl(q,A) as 
a function of A for anisotropy parameter a — 2.5 for L = 0, 
L — 1, and L = 2 (scaled down by a factor of 10). Bottom: 
The total energy difference AE2(a, A) as a function of A for 
various values of the Fermi sea anisotropy parameter a for 
the second excited LL. 



We are interested in the potential and total energy dif- 
ference between the isotropic (a — 1) and the anisotropic 
case (a > 1), namely, AVz,(a, A) = Vl(1,A) — Vt(a, A), 
and AEl(oi, A) = El(1,\) — Ei(a,\), as a function of 
A and the anisotropy parameter a and for L = 0, L = 1 
and L — 2. In the interesting region, the energy dif- 
ference between the nematic and the isotropic state can 
be small relative to the energy scale e 2 /lo, so high ac- 
curacy may be required. The effective potential for high 
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LL and for small values of A changes rapidly at small 
distances and oscillates at large distances. We have used 
an adaptive mesh to incorporate these multi-scale oscil- 
lations accurately for up to L — 2. For higher Landau 
levels it becomes increasingly more difficult to carry out 
an accurate calculation due to the fact that these oscil- 
lations become increasingly more rapid. 

Our calculated pair distribution function for the 
isotropic accurately reproduces the pair distribution 
function reported in Ref. I22I . The potential energy dif- 
ference AVz,(a,A) for various values of a is calculated 
and is shown in Fig. QJtop), in units of e 2 /(elo) (where 
Iq = ^Jhc/eB). Notice that for the case of the LLL and 
for the first excited LL, AVl(oi, A) < 0, i.e., the isotropic 
state is energetically favorable for all values of A and a. 
However, as it is illustrated in Fig. [IJtop), for the case of 
the second excited LL, AV%{a, A) > 0, for all values of a 
and for some range of the parameter A, the anisotropic 
state can be energetically favorable provided that the en- 
ergy loss due to the anisotropy of the Fermi surface is not 
larger than the potential energy gain. 
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FIG. 2: The total energy as function of the anisotropy pa- 
rameter a for various values of A. 



In the single-LL approximation^ the kinetic energy 
of the isotropic state is quenched. We can estimate the 
kinetic energy difference between the isotropic and the 
anisotropic case by ignoring the Landau level projection 
operator and, thus, writing the wave function as \& = 
where $ is the non-interacting Slater determinant and F 
the Jastrow part. The kinetic energy contains terms in 
which the operator (V — A) 2 acts on F. This term gives 
the same contribution of Huj c /2 in both isotropic and 
anisotropic case. Therefore, the main difference, com- 
ing from the term |F , | 2 $*V 2 $, is due to the difference in 
shape of the Fermi sea. This leads to the following ki- 
netic energy difference between isotropic and anisotropic 
Fermi sea: AK ~ — h . fc f . In Fig. QTbottom) we 

4m* 2a G 1 — ' 

present the total energy difference AE^a, A) between 
the anisotropic and the isotropic state. Notice that the 
nematic state is energetically favorable relative to the 
isotropic below A c ~ 0.6. 
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FIG. 3: The optimum values of the total energy of the ne- 
matic and isotropic state as a function A is compared with 
the results of the Hartree-Fock approximation obtained for 
the same values of A. 



Hartree-Fock energies of the stripe and bubble charge- 
density-wave energy have already been reported j 23 i 24 
however, they are only available for the case of A = 0. 
In order to compare the energy of the nematic state with 
that of the ordered stripe and the bubble states for a 
non-zero value of A, we carried out a Hartree-Fock calcu- 
lation using the method outlined in Refs. I23U24I . In Fig. [5] 
the total energy (apart from a common constant value of 
hu> c /2) is compared with the results of our Hartree-Fock 
calculation for finite values of A. In addition, in Fig. [3J 
the minimum energy with respect to the anisotropy pa- 
rameter a is compared to the minimum Hartree-Fock en- 
ergy value with respect to the uniaxial anisotropy pa- 
rameter e (the lattice constants of the uniaxial Wigner 
crystal are given in terms of e as a\ = V3a/2\/l — e and 
(i2 = VI — ea/2). Notice that there is a critical of A, 
namely, A c ~ 0.4 below which the nematic phase is ener- 
getically favored. 

A value for the parameter A can be estimated using 
the calculation presented in Ref. [I?]. Using the value 
of the 2D electron density for these materials^ we find 
that A ~ 62A. This corresponds to a value of A ~ 0.34 
in units of the magnetic length (^o = 18lA for the 
value of B ~ 2T which corresponds to the second ex- 
cited LL in the experiments of Refs. 2.3). This value of 
A is less than the critical value A c ~ 0.4 below which 
the nematic phase is energetically favorable as compared 
to the stripe ordered phase (see Fig. [3]). Therefore, we 
conclude that our calculation suggests that for the case 
of the 2DEG in the heterojunctions used in Refs. 0@ the 
quantum nematic state&^ii 1 - may be energetically lower 
than the stripe ordered or bubble phases. Furthermore, 
it is interesting to probe this transition from the nematic 
to stripe to isotropic by either experimentally altering 
the value of A or indirectly by means of an in-plane field. 

So far, because the Hartree Fock calculations predict 
a stripe state, the observed anisotropy in transport was 
taken as a signature of a stripe ordered state. This state 
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FIG. 4: The potential energy of the nematic and isotropic 
states as a function A is compared with the results of the 
Hartree-Fock approximation. 



breaks both translational invariance in one direction and 
rotational invariance. The results of Haldane and Rezayi 
and Yang2£ are also usually interpreted as a stripe state. 
However, the systems which can be done with such an ap- 
proach are very small, and, in addition, toroidal bound- 
ary conditions were used, which break rotational invari- 
ance. Therefore, because of these limitations it cannot 
be discerned if the true ground state is a stripe or a ne- 
matic. In our calculation we find that the optimum ne- 
matic phase corresponds to an anisotropy of a ~ 10 near 
the physically realized value of the parameter A. This 
implies that a nematic state with such large anisotropy 
cannot be distinguished from the stripe state in systems 



with only twelve electrons^ 

There are two possible sources of systematic error in 
the present calculation. First, the use of the FHNC/0 ap- 
proximation to evaluate the distribution function, where 
the contribution of the elementary diagrams is neglected. 
This approximation works very well in low-density sys- 
tems, i.e., where the average interparticlc distance is large 
compared to a hard core diameter. In the present prob- 
lem such a condition is not clearly fulfilled as there is 
only a soft core of size A. The second source of error is 
the fact that we have neglected the projection operator 
and assumed that the unprojected wave function given 
by Eq. [T] is a good approximation to the lowest LL. In 
order to address these concerns, in Fig. |4] we compare 
only the potential energy of the nematic and isotropic 
states with the results of the Hartree-Fock approxima- 
tion obtained for the same values of A. Notice that for 
values of A > 0.6 the results of the FHNC and the HF 
calculation are almost identical. Moreover, the results of 
the FHNC calculation for the isotropic state agree very 
well with those of the HF calculation for all values of 
A > 0.3. This is an indication that the energy difference 
between the nematic phase and that of the isotropic state 
and the stripe state below A ~ 0.5 may not be an arti- 
fact of the difference in the treatment of the two states 
(i.e., the difference between HF and FHNC approxima- 
tions) but rather due to the fact that the nematic state 
for long-range interactions is energetically favorable for 
at least the second excited LL. 
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